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Abstract 

The absorption of infrared and visible radiation from stellar emission spectra by the negatively 
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charged hydrogen ions H is considered. We derive the explicit formulas which can be used to 
determine the total absorption coefficient (per unit volume) for the negatively charged hydrogen 



ions 1 H (protium) and 2 H (deuterium or D ). The computed bound-free and free-free absorption 



coefficients a u and k v can be used to evaluate the actual absorption of infrared and visible radiation 



by the H ion in photospheres of many cold stars with surface temperatures T s < 8,250 K. 
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I. INTRODUCTION 

As is well known (see, e.g., []| - j^]) the negatively charged hydrogen ions 1 H~ (protium) 
and 2 H~ (deuterium or D~) are of great interest in Stellar Astrophysics. Indeed, the nega- 
tively charged hydrogen ions substantially determine the absorption of infrared and visible 
radiation in photoshperes of all stars with temperatures bounded between T max m 8,250 
K (late A-stars) and T min w 2,750 K (early M-stars). The main contribution to the light 
absorption by the negatively charged hydrogen ion H~ comes from its photodetachment 

H~ + hv = H + e~ (1) 

where the notation hv designates the incident light quantum, while e~ means the free elec- 
tron, or photo-electron, for short. Here and below h is the Planck constant. Note that each 
of the negatively charged hydrogen ions, i.e. the 1 H~, 2 H _ (or D - ), 3 H~ (or T~) and °°H _ 
has only one bound (ground) l 1 ^— state with L — 0, where L is the angular momentum of 
this two-electron ion (see, e.g., [4j). The ground l 1 ^— state is a singlet state, i.e., the total 
electron spin S in this state equals zero, and therefore, 2S + 1 = 1. Here and everywhere 
below the notation ^H/^H" stands for the model atom/ion with the infinitely heavy nu- 
cleus. The total energy of the ground l 1 ^— state in the °°H _ ion is -0.527 751 016 544 377 



196 566 8 a.u. [4|, while the total energy of the hydrogen atom with the infinitely heavy 
nucleus (°°H) is -0.5 a.u. Therefore, the binding energy of the °°H - ion is Xi ~ 0.027 751 
a.u. 0.755 143 9 eV and the corresponding ionization frequency z/ x = ^- is located in the 
far infrared part of the spectrum. Note that the energies of 'visible' light quanta are ^2-3 
eV, while the total energy of the ground state of hydrogen atom is 13.65 eV. Light quanta 
with such energies (~ 13.65 eV) correspond to the vacuum ultraviolet region. The total 
energies of the ground l 1 ^— states of other H~-like ions and H-like atoms can be found in 
Table I. 

The final atomic state which arises after photodetachment of the negatively charged 
hydrogen ion, Eq.flT]), includes the neutral hydrogen atom H and free electron e~ . The kinetic 
energy of this free electron can, in principle, be arbitrary. Therefore, it is easy to predict 
that all radiation quanta with the frequencies v > v\ = ^ can produce photodetachment of 
the negatively charged hydrogen ion H~. In other words, all light quanta with v > v\ can 
be absorbed during such a process. This also includes the quanta of visible light. 
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The absorption of radiation by the negatively charged hydrogen ions has been investigated 
in numerous earlier studies (see, e.g., j|3] - [lo| and references therein). In general, the light 
absorption by an arbitrary atomic system is determined by separate contributions from the 
corresponding bound-bound, bound-free and free-free transitions in this system. As men- 
tioned above the negatively charged hydrogen ion has only one bound (ground) l 1 ^— state, 
i.e. there is no contribution from the bound-bound transitions for the H~ ion. The bound- 
free transitions in the H~ ion correspond to the photodetachment, Eq.flT]), while the free-free 
(electron) transitions formally represent the radiation-stimulated electron scattering in the 
field of the neutral hydrogen atom (or inverse bremsstrahlung, for short) 

e~ + hv + H = H + e" (2) 

where the kinetic energy of the final electron differs from its incident energy due to the 
photon's absorption. If the kinetic energy of the both incident photon and electron in Eq.(J2]) 
are small, then the final state of the hydrogen atom coincides with its original (ground) state. 
Indeed, the lowest excitation energy of the ground state of the hydrogen atom is relatively 
high E pa 0.375 a.u. ~ 10.2043 eV. Therefore, we can ignore all possible excitations of 
the central hydrogen atom during radiation-stimulated electron scattering in the field of 
the neutral hydrogen atom. Such an approximation has a good numerical accuracy, if the 
kinetic (or thermal) energies of electrons E e ~ kT are significantly smaller than the lowest 
excitation energy of the ground state of the hydrogen atom (« 0.375 a.u. ~ 10.2043 eV). It 
can be shown that such an assumption is always obeyed in all stars in which the absorption 
of radiation by the negatively charged hydrogen ions is important, i.e. it is true for the 
late A-stars and all F, G and K stars. Note that in the early F-stars and late A stars the 
absorption of radiation by the negatively charged hydrogen ions is consequently outweighed 
by the neutral hydrogen absorption. In the hot O- and B-stars, with surface temperatures 
T s > 28,500 K and T s > 11,000 K, respectively, the absorption of radiation quanta is 
mainly related with the transitions between different bound states in helium and hydrogen 
atoms. The negatively charged H~ ion does not exist at such temperatures as a stable, i.e. 
bound, system. 

On the other hand, the H~ ions play practically no role for the cold M and N stars with 
T s < 2,700 K where the absorption of radiation by various molecular species and atoms 
of metals becomes important. However, in the late F, G and early K stars the absorption 
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of infrared and visible radiation by the negatively charged hydrogen ions is maximal. This 
includes our Sun which is a star of spectral type G2. The great role of the H~ ions for our 
Sun was suggested by R. Wild in 1939 (see discussions and references in {^J and {3J]). 

In this study we analyze the bound-free and free-free transitions in the negatively charged 
hydrogen ions. Our main goal is to improve earlier evaluations of the photodetachment cross- 
section and the total cross-section of the free-free transitions in the H~ ion. In particular, 
we want to re-consider the validity of some approximations made in earlier works where the 
absorption of radiation by the negatively charged hydrogen ions has been evaluated j^] - [lo| . 
Below, we pay a significant attention to the formulas which are used in numerical calculations 
of the absorption coefficient of the H~ ion. Results of different numerical calculations are 
discussed here very briefly, since we want to include such a discussion in our next paper. 

II. ABSORPTION COEFFICIENT OF THE NEGATIVELY CHARGED HYDRO- 
GEN ION 

In general, the radiation absorption coefficient a u of the negatively charged hydrogen ion 
defined per unit volume is written in the form (see, e.g., [lj]) 



where n~ and nn are the spatial densities of the H~ ions and neutral hydrogen atoms H, 
respectively. Also, in this formula p e = n e kT is the electron pressure which corresponds 
to the temperature T = T e and electron density n e . The factor 1 — exp(^—j^^J in Eq.flH]) 
represents the 'negative absorption' of radiation. The spatial densities nn,n~ and n e are 
related to each other by the following equation 



where g\ = 2 and g_ = 1 are the statistical weights of the ground states of the neutral 
hydrogen atom H and hydrogen negatively charged ion H~, respectively. The notation xi 
stands for the ionization potential of the H~ ion. In this equation and in all formulas below 
the notations m e , h and e mean the electron mass, Planck constant and electron's electric 
charge, respectively. 




(3) 




(4) 
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By using Eq. (Pi]) we can express the spatial density of the H ions (i.e. the n factor in 
Eq.Q) and reduce Eq.Q to the form 



a u = n H Pe 



A(2Trm e )l(kT)-2 ^kT' 
This formula is used for all current calculations of the absorption coefficient a v of the 
negatively charged ion H~. In real calculations of a u it is always assumed that p e = 1 
and Tin = 1, i.e. such calculations are performed for the unit electron pressure and per one 
hydrogen atom in unit volume. In this study the formula, Eq.flH]), is also applied for numerical 
calculations of the absorption coefficient a v of the H~ ion. As it follows from Eq. §5§ in order 
to determine the coefficient a u one needs to evaluate the k v and a v coefficients. The third 
unknown value, which can be found in Eq.([5]), is the ionization potential xi of the H~ ion(s). 
In reality, the ionization potentials Xi °f a U H~ ions are now known to very good numerical 
accuracy (see, e.g., [4j and Table I below). Therefore, below we restrict ourselves to the 
analysis of computational formulas used for the k v and a v coefficients. Note that the first 
numerical calculations of the absorption coefficient a v with the use of Eq. (jSJ) were performed 
in jo]. It is also clear that if the total absorption coefficient a u as the function of r in stellar 
photoshpere is known, then one finds the correspponding optical depth at the frequency v 

r+oo 

r v {r) = / a u {y)dy (6) 



where r is the actual geometrical distance measured from the top of the Solar/stellar pho- 
tosphere (oo). 

Now, let us discuss the relations between the absoption coefficients a u and k v from Eqs.flS]) 
and (jSJ) and the corresponding cross-sections used below. Without loss of generality consider 
such a relation for the bound-free coefficient a v and photodetachment cross-section o v of the 
H~ ion. The general relation produces the following formula a v = a v ■ 1 ■ N = a u ^-N^, 
where Na is the Avogadro number, Mh is the mass of one mole of hydrogen atoms and p is 
the density of the stellar (hydrogenic) photosphere. 



III. PHOTODETACHMENT OF THE NEGATIVELY CHARGED HYDROGEN 
ION 

In this Section we derive and discuss the explicit formulas which are used for numerical 
computation of the photodetachment cross-section of the H~ ion. The photodetachment 
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corresponds to the bound-free (optical) transitions, since the final state of one of the two 
electrons is an unbound state, or, in other words, the state from the unbound spectra of the 
H ion. For the two-electron H~ ion it is possible to obtain the closed analytical formula for 
the photodetachment cross- sect ion, if the non-relativistic (or dipole) approximation is used. 
The derivation of this formula is the main goal of this Section. 

Let us assume that the incident (i.e. second) electron was bound to the neutral hydrogen 
atom with the binding energy = — /, where / = %i is the ionization potential of the H~ 
ion. The incident photon has the momentum k and the frequency u (or energy fiuj = uS). In 
the final state the photo-electron moves as a free particle with the momentum p and energy 
6f. Since p is a continuous variable, the photodetachment cross-section is written in the 
form 

da = 2tt | Mj_>/ | 2 8{-I +u-e f )-0- 3 (7) 

where the wave function of the final state is normalized per one particle in unit volume 
V = 1. All wave functions and expressions used in computations of the transition amplitude 
Mj_>j are assumed to have a proper permutation symmetry upon spin-spatial coordinates of 
the two electrons 1 and 2. This allows one to write the formula, Eq. (!7|) , and all formulas 
below in the one-electron form. 

The delta-function in Eq.flT]) is excluded by integrating over the momenta p of the photo- 
electron. Indeed, by using the formula d 3 p = p 2 d | p | do = e | p | dedo and performing 
the integration over de we can remove this delta-function and obtain the following general 
expression for the cross-section 

da = e 2 ^-^ | Mi-+f | 2 do (8) 

The derivation of the analytical formula for the transition amplitude M^f is drastically 
simplified (see, e.g., [is)]) with the use of transverse gauge, i.e. 

M^f = -eA-j^f = -ey^ J J 2)[(e-ai) exp(ik v r 1 )+(e-a 2 ) exp(2k 2 -r 2 )]^(l, 2)d 3 r 1 dh 2 

(9) 

In this equation and everywhere below we shall use the so-called relativistic units where 
h = 1 and c = 1. 

In this study all electrons are considered as non-relativistic particles, while the energy 
of the incident photon hu)(= uS) is assumed to be larger than /, but it is substantially less 



than the energy of electron at rest, i.e. hu <C m e c 2 . This means that the final velocity of 
the photoelectron is small, i.e. we are dealing with the non-relativistic problem. Therefore, 
the Dirac's a matrixes can be replaced by the corresponding velocities Vj = — ^-Vj, where 
j is the electron index. Moreover, in the non-relativistic approximation we can replace the 
factors exp(«kj • r_j) (j = 1, 2) by uinty. Now, it is easy to obtain the following formula for 
the photodetachment cross-section 

d d 



e 2 m e I p 

da,, = 



2huj 



■{J |**.(n,r 2 )(— + — )v& H -(r 1 ,r 2 )d 3 r 1 rf 3 r 2 }| 2 do (10) 



where e is the vector which represents the polarization of the incident photon, \I/^(ri, r 2 ) is 
the wave function of the final state (after photodetachment), while \I/H-( r i> r 2) is the wave 
function of the ground state of the hydrogen ion. In this study the ground state of the 
H ion is approximated by the following (exponential) variational expansion: 



N N 



*h- = *H-(ri,r 2 ) = J2Cifi = X^ ex P( _a * r i _ Pi r 2 ~ 7^12) (11) 

i=i i=i 

where r 12 =| r x — r 2 | is the electron-electron (or correlation) coordinate and Cj (i = 1, 
. . . , N) are the linear variational coefficients of the approximate (or trial) wave function. The 
3 N— parameters ati, 7* (i — 1, . . . , N) are the non-linear parameters (or varied parameters) 
of the trial wave function, Eq.f llip . As follows from the results of our earlier studies (see, 



e.g., (4], the exponential variational expansion, Eq.f llip . provides very high numerical 
accuracy for any bound state in arbitrary Coulomb three-body systems, i.e. for systems 
with arbitrary particle masses and electric charges. For the H~ ions such an accuracy is 
much higher than accuracy which any other variational expansion may provide, if the same 
number of basis functions iV is used. The explicit form of the \E'^(ri,r 2 ) wave function is 



discussed below. Note that from Eq.flTTJ one finds 

r) N 

+ «r) * H - = ( e • r) £ c ^ h + ( e ' r) S CM (12) 



d d r - N ~ - N 

h — )* H - = (e- — . , . , 

OTi dr 2 / v r x > ^ v r 2 - i=1 



After photodetachment of the H~ ion, Eq.fJTJ), in the final state we have the neutral 
hydrogen atom in one of its bound state and free electron which moves as a free particle 
with the kinetic energy E e = §^-, where % = ^ is the reduced Planck constant and m e is 
the electron mass. The corresponding wave function of the final state is represented as a 



product of the bound state wave function of the hydrogen atom and the wave function of the 
free electron. The wave function of the bound (n£m)— state of the hydrogen atom is written 
in the form (see, e.g., |13|) ® n e m (r, 6, 0) = R n£ (Q , r)Y im (Q , 0), where Y £m (0, (p) = Y im (n) is 
a spherical harmonic and R n g(Q,r) is the one-electron radial function. The radial function 
Rni(Q, r ) takes the form (see, e.g., [131 ]) 

* / 



Rne(Q,r) = — 
rn 



\ 



Q(n-£- l)\v2Qr 



(n + 



n 



e+i 



n-i-l 

{£ 

k=0 



-r 



n + . 
2£ + k 



2Qr 



n 



1 



exp 



where Q (= 1) is the nuclear charge, while n and i are the quantum numbers of this bound 
state. Note that this radial function has a unit norm. 

The wave function of the final (i.e. free) electron is represented as the sum of the plane 
wave 4>(r) = exp(zk-r) and converging spherical wave which is usually designated as (r). 
In the non-relativistic dipole approximation used here one finds that the following selection 
rule is applied: transition from the incident S(L = 0)— state of the H~ ion can only be 
into the final P(L = 1)— state. This means that in the final state we can detect either 
the H-atom in one of its s(£ = 0)— states and free electron which moves away in the p- 
wave (the sp— channel of the H~ ion photodetachment), or the H-atom in one of its three 
p(£ = 1)— states and free electron which moves away in the s— wave (the ps— channel of the 
H~ ion photodetachment). The wave function of the final state for the sp— channel takes 
the form 



3z 

if> f (l, 2) = R n0 (Q = 1, n) — (n e • n 2 )$ p i(r 2 

zp 



3i 
2p 



n e ■ n 2 )i? n o(ri)$pi(r 2 



(14) 



where n 2 = — is the unit vector of the second electron and n e = — is the direction of the 
momentum of the outgoing photo-electron. The notation R n i=o{Q = 1, ri) = R n o( r i) stands 
for the radial function of the bound state of the hydrogen atom with quantum numbers 
n > 1 and £ = 0. 

The wave function of the final state for the ps— channel can be written in the following 
form 

V>/(1, 2) = Rni=i(Q = 1, ri)(n e ■ n x )— $^ =0 (r 2 ) = ^(n e ■ ni)i? n i(ri)$p (r 2 ) (15) 

where R n i(Q = = -Rni(^i) is the radial function of the bound state of the hydrogen 

atom with quantum numbers n > 2 and £ = 1. 
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As follows from Eqs. fllOp - ( TT5|) the formula, Eg. ({TO]) , is reduced to the form 

I ("-)•> + («'-)•• I" < 16 > 



where $1 and $2 are the two scalar expressions. In stellar photospheres all hydrogen ions 
are systems free-oriented in space. This means that the expression from Eq. (ll6j) must be 
averaged over possible directions of the unit vectors ei and e 2 which describe the polarization 
of the incident photon. In the incident light wave we always have e x _L e 2 _L n, where the unit 
vector n corresponds to the direction of propagation of this light wave, or the corresponding 
light quantum (photon). An arbitrary scalar expression can be averaged over all possible 
directions of the unit vectors ei and e 2 with the use of the following formulas 

ZiZk = ^{Sik ~ riiUk) (17) 

and 

(a-e)(b-e) = ^[(a • b) - (a • n)(b • n)] = I(a x n)(b x n) (18) 

where x x y is the vector product of the two vectors x and y. Note that the final formula 
includes only the unit vector n which describes the direction of photon propagation. If the 
direction of the propagation of the final photon is not known, then the last formula must 
also be averaged over all possible directions of the unit vector n. This leads to the formula 

r 87T 1 , . , 47T 



(a • e)(b ■ e) = -(a x n)(b x n) = — • -(a • b) = y(a • b) (19) 

For the sp— and ps— channels mentioned above the radial functions of the final states of 
the hydrogen atom can be designated as the | nO) and | nl), respectively. Here the notations 
n and £ designate the principal quantum number n and angular number £ of the hydrogen 
atom. To designate the radial hydrogen functions we do not need to mention the magnetic 
number (s) m, where m = —£, —£ + 1, ...,£— 1,£. In reality, the largest contribution to the 
photodetachment of the H~ ion comes from that final state in which the final hydrogen atom 
forms in its ground state. This corresponds to the sp— channel. The final hydrogen atoms 
can also be formed in the excited 2s(£ = 0) — ,3s(£ = 0) — ,4s(£ = 0)— states. Such channels 
also contribute to the photodetachment of the H~ ion, but their contributions are much 
smaller than the contribution of the photodetachment channel in which the final hydrogen 
atom is formed the ground state. Analogous situation can be found for the ps— channel 
where the contribution from the 2p(£ = 1)— state substantially exceeds contributions from 



other np(£ = 1)— states in hydrogen atoms (here n > 3). In general, the formation of the 
final hydrogen atoms in the highly excited s(£ = 0)— states and p(£ = 1)— states with n > 4 
is very unlikely, since the corresponding probabilities are very small. 

It is also clear a priori that the contribution of the sp— channel in the photodetachment 
of the H~ ion is substantially larger than the contribution from the ps— channel. It is 
directly follows from very large excitation energy of the 2p(£ = 1)— states in the hydrogen 
atom AE ~ 10.2043 eV. Another reason for this is almost exact orthogonality of the 
wave functions of the 2p— state of the hydrogen atom H and ground state of the H~ ion. 
In reality, in the cold stars the ps— channel does not contribute into the photodetachment 
of the H~ ion. However, the contribution of this ps— channel rapidly increases with the 
temperature due to Maxwellian (thermal) distribution of the H~ ions upon temperatures in 
stellar photospheres. It follows from here that for stars with higher surface temperatures 
the ps— channel must play a larger role in the photodetachment of the H~ ions. In other 
words, there is a difference between the photodetachment of the H~ ions in photospheres of 
the late A-stars and in photoshperes of earlier K-stars. Note also that in our non-relativistic 
approximation the rotationally excited channels with £ > 2 do not contribute in the total 
photodetachment cross-section. 



IV. COMPUTATION OF THREE-BODY INTEGRALS IN PERIMETRIC COOR- 
DINATES 



As we have shown above the problem of analytical and numerical computation of matrix 
elements arising in Eq. ffTO]) with the trial wave function defined in Eq. ffTTj) is reduced to 
the computation of the following three-body integrals with the Bessel functions jo (kr 2 ) and 
ji(Ar 2 ) (see, e.g., [14J and |15j). In fact, each of these integrals takes the form 



4>nt{ri)j m {kr 2 ) exp(-ari - (3r 2 - ^r 12 )r 1 r 2 r X2 dr x dr 2 dr 12 



(20) 



where m — 0, 1 and r 12 =| r x — r 2 | is the electron-electron coordinate. Note that the scalar 
ri2 variable varies between the following limits | r\ — r 2 \< r\ 2 < r\ + r 2 . In Eq. (l20l the 
notation <^w( r i) stands for the unit norm wave functions of the bound states of the hydrogen 
atom. The integrals in Eq. (120]) are computed analytically in perimetric coordinates ui,u 2 , %, 



where «i = \{r 2 +ri 2 — ri), u 2 = h( r i+ r i2—r 2 ) and U3 = \{ri+r 2 — r\ 2 ). The three perimetric 



10 



coordinates u±,u 2 , u 3 are independent of each other and each of them varies between and 
+00. 

Let us derive the closed analytical expression for the integrals, Eq. (l20l) with m = 
and m — 1. First, we need to obtain the general formula for the closely related auxiliary 
three-body integral F^^a, f3, 7) which is written in the form 

r n ,jfc,i(a, A 7) = / / r^r l 2 r^exp(-an - 0r 2 - ^r 12 )drxdr 2 dri 2 (21) 



where k > 0, 1 > 0, n > and a + /3>0,a + 7>0 and + 7 > (see below). Analytical 
computation of this integral has extensively been explained in a number of our erlier works. 
Corresppondingly, below we restrict ourselves only to a few following remarks. In perimetric 
coordinates this integral, Eq. fpZTj) . takes the form 



roo roo pod r -, 

r„, fc ,z(a,/3,7) = 2 / / / exp -(a + 0)m - (a + 7)^2 - (/? + l)u 3 (u x + u 2 ) n 
Jo Jo Jo 1 J 

(wi + u 3 ) m (u 2 + u 3 ) l du 1 du 2 du 3 (22) 

where we took into account the fact that the Jacobian of transformation from the relative 
to perimetric coordinates equals 2. The integration over three independent perimetric co- 
ordinates Ui (0 < Ui < 00) in Eq. (l22|) is trivial and explicit formula for the T nj k,i(cx, fl, 7) 
integral is written in the form 

rWa,/3, 7 ) = 2-/!m!n! (23) 

l m n f~ih (~~m\ s~im\ 

n ^ (a + /3) m -™i+^+ 1 (a + ryy-h+m+i^ + 7 )n-„ 1+mi +i 

where C™ are the binomial coefficients (= number of combinations from k by m). In many 
application all three integer numbers k, I, m are relatively small (< 3) and this simplifies 
numerical applications of the formula, Eq.(f23l). 

The parameters a,/3 and 7 in Eqs. fT20|) - (123]) can be arbitrary real and/or complex 
numbers. Furthermore, one of these three parameters/numbers can be equal zero, but the 
three principal conditions for these parameters a + /3 > 0,a + 7 > and (5 + 7 > must 
always be obeyed. If these three parameters are complex, then these three conditions are 
written in the form Re(a + 0) > 0, Re(a + 7) > and Re{j3 + 7) > 0. These conditions are 
needed to quarantee the convergence of the corresponding integrals in relative coordinates. 

Now, the explicit formulas for the analogous three-body integrals which contain the Bessel 
functions jo(kr 2 ) and/or j\{kr 2 ) is derived in the following way. First, consider the integral, 
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Eq.()20l). with the jo(x) function which can also be written in the form jo(x) = In this 
case the integral Eq. (1201) is written as the sum of the following two integrals 

I (n,£;k) = — / f f (f>nt{ri) exp[-arx - (f3 - ik)r 2 - ^r X2 \r x r X2 dr x dr 2 dr X2 (24) 



2ik 
1 

Ihk 



4>ni(ri) exp[-ari - (f3 + ik)r 2 - ^ri^r^^dr^r^r 
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where i is the imaginary unit. Consider the case when the final hydrogen atom is formed in 
its ground state (or (lO)-state). The integral Io(n, £; k) = Iq(1, 0; k) can easily be determined 
in perimetric coordinates. Indeed, in this case (p n e(ri) = 0io(^i) = 2exp(— Ar), where A = 1, 
and the last integral is reduced to the form 

J (l, 0; k) = J J J exp[-(a + A)r 1 - (/? - ik)r 2 - -fr^r^^drxd^dr^ (25) 
1 

%k 

2^ foo roo roo 

= ~~k 

2,1 fOO roo roo 



+ 



exp[-(a + A)r x — (/3 + %k)r 2 - -yri 2 ]riri 2 dridr 2 dr 12 

roo roo roo 

/ / / expf-Xi^i - X 2 u 2 - X 3 u 3 ](ui + u 3 )(ui + u 2 )duidu 2 duz 
Jo Jo Jo 

'2i r(x> f<x> roo 

— / / exp[-Xi«i - - Xtu^Mux + u 3 )(ui + u 2 )duidu 2 du 3 
k Jo Jo Jo 



where X 1 = a + 7 + A, X 2 = (3 + 7 - ik, X 3 = a + (3 + A - ik and X* 2 = (3 + 7 + ik, XI 
a + (3 + A + ik. Note that each of the last two integrals is written in the form 



roo roo roo | 2 

J(a,b,c)= / / / exp[— ax — by — cz](x + z)(x + y)dzdydz = — — \— 
Jo Jo Jo ' abc a 



111 
+ - + - + -] (26) 
ab ac be 

Due to some hystoric reasons this integral is often called the Coulomb three-body integral. 

The closed analytical expression for this integral allows one to obtain analytical formulas 
for all integrals, Eq. (|20l) . which include the | n£) radial functions of the final hydrogen atom. 
In reality, such integrals are computed as the partial derivatives from Eq. ( 1261) in respect with 
the parameter a. This solves the problem of analytical/numerical calculations of the three- 
body integrals from the exponential functions of the three relative coordinates, Eq.( l20|) . 
which include the Bessel function jo{kr). 

Analytical formulas for the analogous integrlas with the Bessel function ji(kr) are easily 
derived from Eq.(|24l) - ( |25l) with the use of the relation 

. /, s 1 djpjkr) 

Ji(kr) = ---^T (27) 
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known for these two Bessel functions. Analytical and/or numerical computation of the 
corresponding derivatives and integrals is straightforward. In general, all analytical trans- 
formations of formulas are easily performed with the use of different plat forms for analytical 
computations such as Mathematica and/or Maple (see, e.g., 16[ and [17j). Note also that all 



integrals arising during this procedure are regular, since each derivative in respect to k gives 
an extra power of r 2 which compensate the presence of this variable in the denominator of 
the right hand side of Eq. 



V. A SIMPLE METHOD TO EVALUATE THE PHOTODETACHMENT CROSS- 
SECTION 

The negatively charged hydrogen ion H~ is the bound two-electron atomic ion, i.e. N e = 
2, while the nuclear charge Q equals 1. As is well known (see, e.g., j^] - [lO] and references 
therein) the long range one-electron asymptotic of the atomic wave function is written in 
the form 

| #( r ) |= Cr?' 1 exp(-tr) (28) 

where C is some numerical constant, Z — Q — N e + 1, while t = \f21\ and I\ = \i is the 
first ionization potential which corresponds to the dissociation process = H + e~. For 
the H~ ion one finds in the last equation Q = 1, N e = 2, and therefore, Z = and the long 
range asymptotic of the wave function is represented in the Yukawa- type form | ty(r) \ = 
^exp(— tr), where C is a constant which must provide the best correspondence with the 
highly accurate wave function of the H~ ion at large r. The highly accurate wave function 
of the H~ ion is assumed to be known. The explicit formula which is used to determine the 
constant C is C =| \l/ H -( r ) I rexp(tr), where \l/ H -(r) = ^n-(r 32 = r,r 31 = 0,r 2 i = r). As 
follows from the results of actual computations the factor C does not change for relatively 
large interval of variations of r, e.g., for r bounded between 7.5 and 35 atomic units. The 
determined numerical value of this constant C for the negatively charged hydrogen ions 
varies between w 0.12595 (°°H-) and w 0.12605 ^fr). 

The photodetachment cross-section a(H~) of the H~ ion is written in the form (in atomic 
units) (see, e.g., js]) 



64ir 2 C 2 aa 2 p 3 _ 6An 2 C 2 a p 3 
a( ' ~ 3 {p 2 + 7 2 ) 3 ~ 3 {p 2 + 7 2 ) 3 ' ( ' 
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where ao is the Bohr radius and 7 = \f2I\. In the usual units {cm 2 ) the last formula takes 
the form 



3 

a(H-) = 2.7032847 • 1(T 16 ■ C 2 ■ . - P „. , cm 2 . (30) 

yP ~\~ 7 



This formula with the known value of C can directly be used in calculations of the photode- 
tachment cross-section a(H~). This formula essentially coincides with Eq.(8) from [7j. In 



a(H-) = 4.28840 ■ 10~ 18 ■ - P cm 2 , (31) 



our case, C = 0.12595, and therefore, from the last formula one finds 

{p 2 + 7 2 ) 3 

where p —\ p | is the momentum of the photo-electron (in atomic units) and 7 = 
0.23558869473885 (see Table I). 

The overall accuracy of this simple formula, Eg. (I31[) . is outstanding. This remakable fact 
can be explained by using the following inequality a< A, where a ~ a is the effective radius 
of the atomic H~ ion (ao is the Bohr radius), while A is the wavelength of the photon which 
produces photodetachment of the H~ ion. For the H~ ion we have a ~ ao ~ 0.529A and 
A ~ 1-10 5 A Briefly, we can say that the incident photon can see only an outer asymptotics of 
the three-body wave function. Briefly, we can say that all details of the wave function of the 
H~ ion at small inter-particle distances are not important to describe the photodetachment 
of this ion. Recently, this simple method was often used in various applications related to 
the H~ ion. However, some important details of the photodetachment of the H~ ion cannot 
be investigated with the use of only this simple method. 



VI. INVERSE BREMSSTRAHLUNG IN THE FIELD OF THE HYDROGEN 
ATOM 

In this Section we discuss the electron scattering stimulated by the electromagnetic radia- 
tion in the field of the heavy hydrogen atom, Eq.([2]). In classical physics this process is called 
the inverse bremsstrahlung. In stellar photoshperes of hydrogen stars inverse bremsstrahlung 
is responsible for absorption of infrared radiation with A > 12,000 A. Furthermore, it is very 
important to describe the total absorption of infrared radiation in photospheres of many cold 
stars. 

As mentioned above at small energies of the incident photon hv < 5eV the original 
(ground) state of the hydrogen atom does not change. This means that for photons of small 
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energies we can ignore all possible excitations of the hydrogen atom, i.e. atomic transitions 
from the ground state of the H atom into other excited states. The incident photon is 
absorbed exclusively by the electron which moves in the field of the neutral hydrogen atom. 
In this case the electron scattering at the hydrogen atom (H + e~ + hv = H + e~) can be 
considered as a regular (but non-elastic!) scattering of the electron in the field of the heavy 
hydrogen atom. In the lowest order approximation we can replace the actual hydrogen atom 
by its polarization potential, which is a central and local potential V(r) — — %, where a is 
a positive constant. In particular, for the °°H atom in atomic units one finds a = |. Below, 
we shall consider the electron scattering in the field of the following central potential 

V ^^ b) = -^fw (32) 

which analytically depends upon the two real parameters a and b. We need to determine the 
cross-section of inverse bremsstrahlung for an electron which moves in the field of the central 
potential V(r; a, b), Eq. (j3"2|) . By using different numerical values for the a and b parameters 
we can obtain a variety of such cross-sections. The main interest, however, will be related 
to the cases when a ~ 4.5 and 6 ~ 0, respectively. 

First, let us obtain the incident and final wave functions which are the stationary one- 
electron wave functions of the continuous spectra of the radial Schrodinger equation (see, 

n 

e.g., [13[) of the electron moving in the field of the central potential, Eq. (l32p . The solution of 
one-electron Schrodinger equation can be found by representing the unknown wave function 
$?(r) in the form \t^(r) = Ag(r) exp(i5e(r)). Note that in the field of any central potential the 
angular momentum £ is a good (i.e. conserving) quantum number. It reduces the original 
problem to the solution of the following differential equation 

dS e {r) 2V(r;a,b)r 12 
dr = 7 [3e{kr) cos 5 e (r) - n e (kr) sin d £ (r)j (33) 

where the phase 5g(r) of the wave function is assumed to be explicitly r— dependent and 
8i(r — 0) =0 (initial condition). With this initial condition the last equation is reduced to 
the form of the following integral equation 

2 f + °° r ]2 

$e( r ) = — t / V (r; a, b) ji(kr) cos 8g(r) —ni(kr)sm6e(r) dr (34) 



k JO 

The equations Eq.( l33|) and Eq.( l34"j) are the main equations of the variable phase method. 
The method of variable phase has been developed by V.V. Babikov and others in the middle 
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of 1960's [18j (see also references therein). Each of these two equations Eqs. (l33l - ( 134)) can 
be solved to high accuracy by using various numerical methods. In general, if we know the 
phase function <5^(r) for each r, then it is possible to determine the amplitude Ag(r) of the 
wave function and the wave function \I/fc(r) = Ap(r) exp(i8?(r)) itself. The corresponding 
equation for the amplitude Ag(r) takes the form 

dA e (r) A e (r)V(r;a,b) 



18| 



dr k 



ji(kr) cos 8e(r) — ng(kr) sin 8g(r) je(kr) sin 8g{r) 

+ng(kr) cos 8i(r) (35) 



Note that the variable phase method works very well for central potentials which rapidly 
vanish at r — > +00, e.g., for | V(r) |~ r~ n , where n > 3 at r — > +00. This means that we can 
apply this method to determine the accurate wave functions of the continuous spectra in the 
field of the polarization potential, Eq. (!32|) . created by the neutral H atom which is infinitely 
heavy. It appears that the variable phase method provides a large number of advantages 
in applications to the problems related with the absorption/emission of radiation from the 
states of continuous spectra. 

Note that in the case of inverse bremsstrahlung the energy conservation law is written in 
the form 

h uj = — - , or nuj = — (36) 

where ki and kf are the wave numbers defined as k a =| k a |= and p a is the momentum 
of the particle a and k a is the correspponding wave vector. In Eq. fl36l) the notation h 
is the reduced Planck constrant, while m e is the electrons mass. From Eq. (l36|) one finds 



k f = + H%~ UJ - us i n S the variable phase method described above we can obtain the 
initial and final state wave functions, i.e. ^/^(r) and \E f / C/ (r), respectively. In the variable 
phase method these wave functions are represented as the products of their amplitudes 
and phase parts, i.e. ^ t,kX r ) = A t;kA r ) ex P^8 i]kt (r)) and ^i ]kf (r) = Ag. kf (r) exj)(i8 i;kf (r)). 
These two wave functions are used in the following computations of the cross-section of 
inverse bremsstrahlung. The explicit formula for the cross-section of inverse bremsstrahlung 
is (in relativistic units) 

,2 

(27i) i m e p i 



da ^Pf = ^oJviL ~ I e *P'^/ I 2 do k d 3 p f (37) 



where Pi =\ Pi \ is the 'absolute' momentum of the incident electron, P/ is the momentum 
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of the final electron and Pi_>/ is the following matrix element 

P*-/ = -»/ ^V^cft- (38) 

where ^/^ and \1// C/ are the incident and final one-electron wave functions (see above). 

Note that the wavelength A of the light quantum which produces photodetachment of 
the H~ ion is > 12,000 A, while the effective 'scattering zone' R s for the process, Eq.(j2]) 
is less than 50 A. Therefore, we can apply the approximation based on the fact that the 
dimensionless parameter ^ < 1. This approximation allows one to replace the r— dependent 
phase factor Sg.^r) in the wave function by its asymptotic value ^.^(+00) determined at 
the infinity. In other words, the phase of the wave function does not depend upon r, i.e. it is 
a constant. The wave function ^.^(r) is represented in the form ^/^.(r) = Ag(r) exp(z^), 
where 5e = 5e(r = +00) is a constant (or almost a constant). This drastically simplifies 
numerical computations of the integrals in Eq. (l37"l) and Eq.(l38l). 

VII. CONCLUSION 

We have discussed the absorption of infrared and visible radiation by the negatively 
charged hydrogen ions H~. Such an absorption of radiation plays an extremely important 
role in Solar/Stellar Astrophysics. In general, the absorption of radiation by the H~ ions 
is mainly related with the photodetachment of these ions. The absorption of radiation by 
electrons scattered in the field of the neutral hydrogen atom(s) also contributes to the total 
cross-section at very large wavelength A > 12,000A (infrared radiation). The photode- 
tachment cross-section of the H~ ion has been determined with the use of highly accurate 
variational wave functions constructed recently for this ion. We explicitly discuss all details 
important for calculations of the photodetachment cross-section with the use of highly ac- 
curate (or truly correlated) wave functions of the three-body H~ ion. To evaluate the total 
cross-section of electron scattering in the field of the polarization potential of the neutral H 
atom we apply the variable phase method. The agreement between our cross-sections and 
results obtained in earlier works for different wavelength can be considered as very good. 
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TABLE I: The non-relativistic total energies E and ionization potential xi °f the hydrogen nega- 
tively charged ions (isotopes) in atomic units. 



ion E xi 

X H~ -0.527 751 016 544 377 196 567 0(2) -0.027 751 016 544 377 196 567 
T~ -0.527 649 048 202 999 950 354 2(2) -0.027 739 991 660 945 743 753 5 
D- -0.527 598 324 686 477 854 280 7(2) -0.027 734 509 437 460 880 345 7 
X H- -0.527 445 881 114 178 893 391 3(2) -0.027 718 041 397 713 212 403 8 
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